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ABSTRACT 
We  consider  the  problem  of  finding  a  polygon 
nested  between  two  given  convex  polygons  that  has 
a  minimal  number  of  vertices.  Our  mam  result  is  an 
O(nlogfc)  algorithm  for  solving  the  problem,  where 
n  is  the  total  number  of  vertices  of  the  given 
polygons,  and  it  is  the  number  of  vertices  of  a 
minimal  nested  polygon.  We  also  present  an  0(n) 
sub-optimal  algorithm,  and  a  simple  0(,nJc)  optimal 
algorithm. 

1.    InlroducLion 

We  provide  an  efTicient  algorithm  for  the  follow- 
ing problem:  given  two  convex  polygons  P  and  Q 
such  that  Q  is  contained  in  F.  determine  a 
minimum  vertex  polygon  K  that  contains  <?  and  is 
contained  in  P  A  polygon  K  is  called  nested 
between  P  and  Q  when  it  circumscribes  Q  and  Is 
inscribed  in  P-  The  problem  was  originally  posed  by 
Victor  Klee  lor  polylopes  in  arbitrary  dimensions;  to 
the  authors'  knowledge  this  remains  an  open  prob- 
lem.   For  the  two-dimensional  case,  we  present  an 


algorithm  that  runs  in  0{nJ.ogk)  time  where  n  is  the 
total  number  of  vertices  of  P  and  Q.  and  k  is  the 
number  of  vertices  of  K. 

The  problem  of  determining  a  minimum  vertex 
nested  polygon  belongs  to  the  general  area  of  polyg- 
onal approximations.  This  area  has  been  extensively 
studied  recently  due  to  its  numerous  applications 
in  robotics,  stock  cutting,  collision  avoidance  and 
Computer  Aided  Design  problems  (see  [CY],  [DA], 
[DB],  [KL]).  Unlike  previous  investigations,  our  goal 
is  to  minimize  the  combinatorial  complexity  of  the 
approximation,  rather  than  a  continuous  measure 
such  as  area. 

Our  paper  is  organized  in  five  sections  Section 
2  introduces  some  notation  and  basic  lemmas.  Sec- 
tion 3  provides  a  straightforward  linear  algorithm 
that  yields  a  nested  polygon  wiLh  at  niosl  one  ver- 
tex more  than  the  minimum.  Using  this  algorithm, 
a  simple  0(Tik)  algorithm  is  developed  This  algo- 
rithm is  subsequently  improved  to  run  in  O(nlog*:) 
time,  which  is  the  main  result  of  the  paper  Finally, 
Section  5  concludes  with  some  remarks  and  direc- 
tions for  further  research. 


2.    Dasic  Letninas 

In  this  section  we  present  basic  lemmas  that 
characterize  minimum  vertex  polygons  and  help  us 
discrelize  the  search  routine. 

A  nested  polygon  is  called  TniniTnal  if  it  has  the 
minimum  number  of  vertices  among  all  polygons 
nested  between  P  and  Q  Henceforth,  we  will  use  P 
and  Q  to  represent  the  boundari«  of  the  respective 
polygons.  The  vertices  of  the  polygons  are  assumed 
to  be  indexed  counterclockwise.  Let  A  be  the 
closed  region  (the  annulus)  of  the  plane  bounded  by 
P  and  (?.  . 

Lemnia  1:   Every  minimal  polygon  is  convex. 
Proof:  Assume   to  the  contrary  that  there  exists  a 

minimal  polygon   K-{-u^ i/^)   that  is   not  convex. 

Let  v^  be  a  non-convex  vertex.  Convexity  of  P  and  Q 
assures  that  the  segment  v«-,v<n  lies  completely  in 

A.    Thus   t/o Vj_,,u<n,....t(t.,    is   a    nested    polygon 

with  one  fewer  vertex.  This  contradicts  the 
minimality  of /C.    Q.E.D. 

A  suppoT-iiny  line  segment  is  a  directed  seg- 
ment in  A  that  supports  Q  on  its  left  and  has  both 
its  endpoints  on  P.  For  any  point  a  on  /-"  .  let  I,  be 
the  unique  supporting  segment  ab  ,  directed  from  a 
to  6.  Define  R,  to  be  the  region  //,  nA.  where  H,  is 
the  closed  right  half-plane  determined  by  the  line  i„ 
(Kig.  1). 

Lemma  2:  For  any  aeP,  R,  contams  at  least  one 
vertex  of  any  minimal  polygon 

Proof:  If  /?,  does  not  contain  any  vertex  of  K.  then 
there  must  exist  two  adjacent  vertices  of  K,  i  and 
y,  that  are  joined  through  /?,.  As  both  x  and  y  are 
to  the  left  of  (,,  it  is  impossible  for  the  edge  xy  to 
intersect  /?,.  Q.E.D. 

Define  a  siipporiing  poLygon  (v^.vg..  .  .v^)  as 
one,  all  of  whose  edges,  except  for  perhaps  the  last 


Pig.  1    A  supporting  segment  ab. 

one  (vttij,  are  supporting  segments.  Let  5,  be  tlie 
supporting  polygon  with  vj  =  a.    For  an  edge  U(U(,i 

of   the   supporting   polygon  S,  =(ii,,v, v*),   define 

the  Q-contact  as  the  vertex  of  Q  that  supports 
V(V,n  and  the  P-conlacl  as  the  edge  of  P  on  which 
Vi  lies.  In  case  of  degeneracies,  we  choose  the  ver- 
tex of  Q  and  the  edge  of  P  with  higher  indices. 
Clearly,  a  supporting  polygon  of  fc  edges  has  at  least 
k-l  ^"Contacts  and  k  P-contacts.  Also,  S,  can  be 
completely  and  uniquely  specified  by  aef,  the  start 
vertex,  and  all  its  contacts  with  P  and  9- 

Lemma  3.  For  any  aeP,  S,  has  at  most  one  vertex 
more  than  the  minimum. 

Proof.     5,    =   (v,, .v^)   has.   by  definition,   at  least 

(jfc-1)  supporting  segments,  namely  1^,.  'v,,,  By 
Lemma  2,  each  /?„,,  defined  by  („,,  where  Q<i^(k  -\), 
must  contain  a  vertex  of  the  minimal  polygon. 
Hence  every  minimal  polygon  must  have  at  least 
(fc-1)  vertices.  Q.E.D. 

In  the  following  lemma,  we  establish  that  the 
search  domain  for  the  minimal  polygon  can  be  res- 
tricted to  one  /?„ 

Lemma  4:    Any   minimal   polygon   K=-{vt.    .  .v^)   can 

be  transformed  into  an  5,  for  some  aeP. 

Proof      The  transformation  is  accomplished  in  two 

steps: 


(1)  Force  all  vertices  of  K  to  lie  on  P.  CadJ  this 
intermediate  polygon  K". 

(2)  Identify  one  vertex  of  A'  as  a  and  rotate  all 
the  sides  of  K"  until  all,  except  perhaps  the 
last  one,  become  supporting  segments,  pro- 
ducing the  polygon  K". 

Step  1:  (Fig.  2a)  Let  v,.  v^,  and,  wj  be  three  con- 
secutive vertices  of  K  such  that  vj  does  not  lie  on 
P.  Let  the  directed  segment  Vji/j,  when  extended, 
intersect  P  at  vj  Replace  vj  by  vj  Convexity  of  P 
and  Q  assures  that  v^uj  lies  completely  in  A. 
Minimality  of  K  assures  that  neither  a  colUnearity 
nor  a  non-convexity  can  arise  from  this  alteration. 
This  procedure  can  be  iteratively  applied  to  move 
all  the  vertices  of  A"  onto  P  The  number  of  vertices 
of  the  resulting  polygon  K"  is  clearly  the  same  as 
that  of  K. 


establishing  the  lemma. ^.^'.D. 


fi^ .  2a.     Transfortriatiori  lo  a  supporting  polygon 
(step  1). 


Step  2:  (Fig.  2b)  Without  loss  of  generality,  set 
v,=i/,=a.  Rotate  v^v^  about  uj  bntU  it  becomes 
coincident  with  i,.  Replace  vj  with  the  other  end- 
point  of  la,  uj  Repeat  the  procedure  for  vjt/j  and 
so  on  up  to  ut_|Ki.  Again,  the  minimality  of  K" 
ensures  that  no  collinearities  or  non-convexities 
can  occur  in  this  procedure. 

After  the  completion  of  Step  2,  K"  is  supporting, 


Fig.  2b     Transfortnalion  lo  a  supporting  polygon 
(step  2). 


Lemmas  2,  3  and  4  are  sufficient  lo  establish 
that  a  minimal  polygon  can  be  found  by  examining 
all  5,s  where  i  is  a  point  on  P  in  /?„ ,  for  any  aef. 

3.    Algorilhms  for  finding  S^  and  minimal  polygons 

In  this  section,  we  develop  an  algorithm  to  find 
Sa.  given  a  point  a-  on  P,  and  use  it  lo  find  a 
minimal  polygon. 

ALgonthm  for  consirricting  S^ 

5a  IS  constructed  by  first  identifying  u,  with  a 
and  then  wrapping"  around  Q  with  supporting  seg- 
ments iv^  ,  i„j  .  i„j ^  The  Inlersection  of  l^^  with  P  is 
taken  to  be  i/,,,,  for  l^i^*-!.  Only  Ihc  last  edge 
may  not  be  a  supporting  segment,  m  which  case  it 
IS  forced  to  be  v^v^. 

The  segment  iv,  can  be  found  in  0(n)  steps  by  a 
linear  search  around  Q  for  the  ^contact  q,  and 
around  P  for  vj  Subsequently,  each  new  edge,  ^  is 
found  in  time  linear  in  the  number  of  edges  of  P  on 
the  boundary  of  /?„  and  vertices  of  Q  between  the 
9-contacts  of  i„,  ^  and  fv,,  q,.,  and  ijt  Therefore,  the 
time  complexity  ol  this  algorithm  is  0(n).  This 
establishes  our  claim  that  a  sub-optimal  nested 
polygon,    with    at   most   one   vertex   more    than    the 


minimum,  can  be  conslrucled  in  linear  lime. 

In  the  following  we  describe  how  to  obtain  a 
minimal  polygon  by  modifying  an  arbitrary  sup- 
porting polygon.  Let  a  point  zEP  be  called  a  con- 
tact c/iange  point  if  S,  has  at  least  one  contact 
different  from  Sy,  where  ye/'  is  a  point  in  the 
neighborhood  of  z  immediately  preceding  it  (clock- 
wise). 

Intuitively,   a   minimal   polygon   is  obtained   by 
starting  with  S,  =5.  and  rotating  5,  around  in  A  by 
sliding     z     along     P     and     checking     Lf     the     non- 
supporting  segment  v^v ,  ever  collapses  to  a  point. 
This  search  is  discretized  by  c^olTTputing  5,  only  for 
those    pomts   i    along    P    that   are    contact   change 
points.     Furthermore,    starting    with    an    arbitrary 
supporting  polygon,  S,,  we  need  to  search  through 
only   those   contact  change  points  that  lie  on   the 
outer  boundary  of  /?,.    It  may  be  the  case,  however. 
that   a  minimal   polygon  5,    is  achieved   at  some  x 
that  lies  beliueen  contact  change  points.    To  locate 
such  minimal  polygons,  we  use  a  set  of  projection 
functions  that  relate  the  position  of  any  vertex  v,, 
m   particular  vt,  to   that  of  t/,    These   functions,   as 
described    below,    are   simple    polynomial   quotients 
and   have    Ihcir  coefTicients   dependent  only  on   the 
current  P-  and  ^-contacts.    As  the  contacts  of  S,  do 
not  change  while  z   is  moved  between  two  adjacent 
contact    change    points,     these     functions    do     not 
change    either      Therefore,    two    adjacent    contact 
change    points    also    determine    the    range    for    the 
applicability  of  these  functions,   within  which  they 
can  be  used  to  detect  a  coincidence  of  v^  with  v,. 

ProjecHon  F-nnctions 

Let  5„|  =  ("i.  vj.  —  nt)  be  a  supporting 
polygon.  Let  the  ^-contact  of  v^v^  be  g=(7..g„).  Let  i 
and  e  be  vectors  parallel  to  the  edges  of  P  contain- 


ing V|  and  vt  respectively.  Assume  a  coordinate  sys- 
tem with  its  X-axis  aligned  with  i.  Using  i,  and  zj  as 
parametric  representations,  the  positions  of  v,  and 
vt  can  be  represented  as  follows: 

y,  =  X, I 
»2  -  d  +  x;  e 

where  d  is  some  constant  vector  (see  Fig.  3). 


/'Xg.  3      Computing  projection  functions. 

Let  r  and  s  be  the  points  on  A'-axis  that  coincide 
with  the  perpendicular  projections  of  q  and  vj, 
respectively.  From  the  similarity  of  the  triangles 
Av,qr  and  t\v^v2S,  we  obtain 

1y/  {1t-^\)  -  'z  siniS/  (d-x  i+ijcosd). 
Solving  for  ij  gives 

ij  =  6(d -x,)/ sim3(gi -(j^cDii3-i|) 
which  IS  equivalent  to 

i2  =  (c,  +  C2i,)/((;3  ■(-  r;,!,)  (I) 

for   constants   c ,.  Cj.  Cj  and  c,    that  depend    only   on 

the  contacts  of  v^v^-    In  the  same  way,  Ihe  position 

of  uj  is  functionally  related  to  that  of  vj.  i.e  . 

I3  =  (d,  +  d2X2)/(d3  +  djij)  (2) 

for   constants   d,.  dj.  d,  and  d,    tliat   depend    only   on 

the  contacts  of  i/jvj.  Substitution  of  ij  from  (1) 
into  (2)  gives,  after  simplification,  a  composed  rela- 
tionship between  13  and  z,: 

z,  =  (a,  +  gjZ|)/(e3  1- «,z,)  (3) 

where    a,.  Bj.  g,  and  b,    are    constants    that    depend 


only  on  the  contacts  between  vi  and  V3.  In  the 
same  way  i,,  i,,,  ...  it  can  also  be  expressed  as 
functions  of  z,.  We  will  write  these  functions  as  Xj  = 
/<(ii).  As  the  degree  of  /,  is  independent  of  i, 
these  functions  can  be  evaluated  in  constant  time 

We  now  describe  how  to  compute  the  next  con- 
tact change  point  for  a  given  supporting  polygon. 
ComputaiioTi  of  Next  Contact  Change  Point 

Given  a  supporting  polygon  5v,  =  ("i.  ^i vt), 

we   need    to   determine    the   edge   of   5„|    that   first 

changes  its  contact  if  v,  is  slid  along  P  while  the 
polygon  is  constrained  to  remain  supporting.  The 
corresponding  position  of  t/,  is  called  the  next  con- 
tact change  point  of  5„  .  We  flrst  compute  the  edge 
that  undergoes  this  contact  change  and  then  use 
the  projection  functions  to  compute  the  next  con- 
tact change  point. 

For  an  edge  v^Vi^,  of  S^^.  let  iji  be  its  ^contact. 
Rotate  V(ii,,i  (extended  to  P  on  both  sides)  about  71 
until  either  the  edge  v^v^^^  becomes  flush  with 
7i7<n.  or  either  v,  or  u,,,  coincides  with  a  vertex  of 
P  Let  the  corresponding  position  of  v,  be  denoted 
by  vf  (see  Fig.  4). 

if  V,  is  moved  from  its  current  position  to  v' 
while  the  polygonal  chain  of  vertices  between  w,  and 
v,n  is  constrained  to  remain  supporting,  at  least 
one  edge,  v^vn,,  changes  its  contacts.  Let  1//  be  the 
first  position  between  v,  and  vf  such  that  moving  v, 
to  V('  forces  some  edge  among  v^v^.  v^v^,  ....  ,  vji/j^, 
to  change  contact.  Our  goal  is  to  compute  vi-,  and 
the  edge  that  changes  its  contact  first.  Assume 
inductively  that  we  l<now  tj,*_,  Then  u,'  is  the  more 
clockwise  of  the  two  points  v{  and  the  intersection 
of  the  extension  of  v^_^q^  with  P-^KSg  4).  This  calcu- 
lation  13  performed  by  the  for  loop  in  the  following 
pseudo-code: 


'i+l 


Pig.  4      Cortiputing  next  contact  change  point. 

ALgorithm  for  computing  a  minimal  polygon 

1.  Choose  a  point  aep.  (Let  a  be  a  vertex  of  P.) 

2.  Compute  the  initial  polygon,  S,  and  set  u,  =  a 

3.  Compute  f/d,  the  set  of  projection  functions. 

4.  whiie   u,e/?3  do  begin 

iCompute  next  contact  to  change  j 

5.  Compute  v^;        i*    .-    2;        vj*-    intersection 
of  v^g,  extended  to  P; 

6.  for  i  =  2  toifc-l     do  begin 

7.  Compute  I/'; 

9  if  Vi  IS  clockwise  of  v,'   then  begin 

9  V('  <-  v{; 

10.  i'  ,.  I 

cndfifi 
11-  v/n  «-intersection  of  1/,'ij,  extended  to  P; 

end  fforj 

12.  Compute  next  contact  change  point  from  i* 
A-': 

13.  Move  f,  to  next  contact  c/iange  point; 

14.  Check  for  the  overlap  of  v^  with  v,  between 
the  contact  changes  by  solving  xt=/t(i|): 

15.  Recompute/,',/,',,, /»' 


end  jwhiiej 

CompLexUy  of  the  Algoriihrn 

Steps  1  and  2  require  0(n).  Slep  3  Lakes  0(nk).  The 
while  loop  of  step  4  is  executed  as  many  limes  as 
Ihe  number  of  contact  changes  that  can  occur  dur- 
ing the  movement  of  v,  throughout  /?..  That  the 
total  number  of  contacts  that  can  change  during 
the  movement  of  v,  throughout  R^  is  0(n)  can  be 
established  as  follows.  Let  the  initial  supporting 
polygon  5,  define  a  set  of  disjoint  regions 
/?,,/?2,...,y?t_i,  one  for  each  supporting  edge  of  5.,  as 
defined  before.  As  v,  is  moved  through  /?i.  the  i"* 
edge  of  the  corresponding  supporting  polygon  will 
have  the  following  as  contacts  throughout  this 
movement: 

(1)  the  edges  of  P  that  lie  iif"SiT>nd 

(2)  the   vertices   of   Q   that   lie   between   the    Q- 
contacts  of  the  segments  defining  /?,  and  /?(n. 

Hence,  all  the  edges  of  S^^  together  can  have  0(n) 
contact  changes  as  v,  moves  through  out  /?,.  Each 
execution  of  the  while  loop  takes  0(k)  time  to  com- 
pute the  ncxl  contact  change,  and  0{k)  time  to 
update  functions  affected  by  this  contact  change. 
Checking  for  overlap  between  v^  and  v,  can  be  done 
in  0(1)  time.  Hence  the  overall  time  complexity  of 
the  algorithm  is  0(nk). 

-1.    O(nlogit)  Algorithm  for  mlninial  polygons 

The  O(njt)  algorithm  for  minimal  polygons  cam 
be  improved  to  0(TUogk)  by  organizing  the  edges  of 
the  nested  polygon  and  the  projection  functions 
hierarchically  into  a  complete  binary  tree  so  that 
the  movement  between  two  contact  changes  can  be 
accomplished  in  0{logk)  rather  than  0(k)  time  We 
shall  assume  throughout  the  following  discussion 
that  k-l^Z*  for  some  integer  d.  In  Section  3.  we 
introduced   the   projection  functions  /,   for  a  sup- 


porting polygon  (vi.vj vj),  which  map  the  posi- 
tion of  a  vertex  t;,.  expressed  as  a  distance  along 
vj's  contact  edge,  onto  a  position  along  vertex  v,'s 
contact  edge.  Now  we  need  the  ability  to  map 
between  any  two  vertex  coordinate  systems  Let 
/,j  be  the  function  that  maps  the  position  of  ver- 
tex V,  to  a  position  along  v^'s  contact  edge,  where 
i>j.  These  functions  are  similar  in  form  and  proper- 
lies  to  /t,  and  in  fact  fi,\=fi. 

Our  main  data  structure  is  a  complete  binary 
tree  with  k~l  leaves  that  are  indexed  by  edges  of 
the  nested  polygon.  The  indices  are  in  increasing 
order  from  right  to  left,   that  is.   leaf  i   represents 

the   edge  fjv,,,   of   the    nested    polygon   (V|.uj v^). 

For  any  node  s  in  the  tree,  let  the  subtree  for  which 
it  is  the  root  be  denoted  by  r, .  A  node  s  In  the  tree 
stores  the  foUowing  information: 

(1)  i(s)  =  irain(r,).  index  of  the  leader,  that  is, 
the  lowest  indexed  (rightmost)  leaf  in  T,, 

(2)  /i.i(.)  where  t=  H-imax(r,).  with  imax(r,)  the 
largest  index  of  a  leaf  of  T,,  [When  i  =  *.  the 
function  /t,i(,)  is  never  actually  used,  and  we 
assume  that  /ti-i  is  the  identity  ) 

(3)  i(s),  index  of  the  first  contact  to  change 
among  the  leaves  of  7',  when  u,(,)  is  moved  coun- 
terclockwise; 

(4)  i(s),  the  distance  by  winch  u,,,,  must  move 
before  i(s)  changes  its  contact. 

(5)  i(s).  a  correction  factor  Ihat  records  the 
distance  by  which  the  leaves  of  T,  have  moved 
when  a  lower  indexed  edge  changes  contact. 

This  notation  is  illustrated  in  Kig    6. 

Let  this  tree  be  called  a  function  tree  and  its 
root  be  denoted  by  r  Let  sj,  and  sr.  respectively, 
be  the  left  and  right  children  of  a  node  v  Given  a 
supporting  polygon  (fi-^s,  ■'"t)  and  the 
corresponding    function    tree,    the    next   contact   to 


change  can  be  computed  in  0(1)  lime:  il  is  deLer- 
mined  by  comparing  the  winners"  from  the  two  sub- 
trees, r,^  and  Trg  within  the  coordinate  system  of 
l(r)  =  i,  i.e..  by  computing 

mm  \fi.,{6(r^))  .  6(r,,)\ 
where  i  =  l(r^)  =  {k-\)/2+l.  Note  that  l{r^)=l  so  that 

both  distances  are  expressed  in  the  same  coordi- 
nate system. 

Hence,  our  main  tasl<  is  to  construct  this  tree 
in  Oinlogk)  time  and  update  it  appropriately  after 
each  contact  change  in  O(logk)  time 

Initial  ConstnLclion  of  Ihe  Function  Tree 

The  tree  is  computed  recursively,  bottom-up 
from  the  leaves.  The  information  associated  with  a 
leaf  node  >  is  easily  computed  in  constant  Lime  as 
follows: 

(1) '(;■)=;; 

(2)  fjti.)  is  computed  as  discussed  in  Section  3: 

(3)  10)=;; 

(4)  i(j)  =    distance  from  Vf  to  i//  (see  section  3 
for  iy,'); 

(5)  A(;)  =  0. 

The  information  associated  with  an  internal 
node  s  is.  then,  computed  from  its  two  children  as 
follows: 

(l)i{s)  =  i{s^): 

(2)  /(.!(.)  =  /i.i(.i)  /i(.j).((.,)  ,  where  i=iTnas{T,): 

(3)  i(s)    is    either    i(st)    or    i(sr)    depending    on 
which  achieves  6{s)  in  (4)  below; 

(4)  (5(s)  =  niin(/,|((5(st)).i5(.i„)),  where  i  =  i(si)  and 
;=^ss). 

(5)  A(s)  =  0. 

Since  the  computation  at  the  leaf  level  can  be  per- 
formed in  0(*)  time,  the  entire  tree  can  be  con- 
structed in  O(kLogk). 


Updating  the  tree  after  a  contact  change 

We  show  that  only  0(logk)  nodes  in  the  tree 
need  to  be  updated  after  a  contact  change  has  been 
performed.  Let  the  last  contact  to  change  be 
indexed  u.  The  leaf  information  associated  with  u 
can  be  recomputed  in  0(1)  time  following  the  ini- 
tialization step  discussed  previously.  Since  the 
relative  ordering  of  the  contacts  in  a  subtree  7",  is 
not  aflected  by  the  change  of  a  contact  that  is  not  a 
leaf  in  7",,  correction  of  the  information  stored  at  s 
would  implicitly  update  all  of  T,-  The  tree  is 
updated  by  recomputing  the  associated  function 
and  adjusting  the  correction  factor  A  of  all  the 
sibling  nodes  along  the  path  from  u  to  the  root  r. 
There  can  be  O(logk)  such  nodes  We  now  show  that 
each  node  can  be  updated  in  0(1)  Lime.  The  correc- 
tion of  f\ji{,)  and  t(s)  can  be  done  by  the  same  com- 
putation as  in  the  initial  construction  of  the  tree. 
The  remaining  updates  are  done  as  follows: 

i5(s)  =  mini/(.y((5(si)  +  a(st)).  i5(sR)  +  6(s;,)i,        where 

i=l{si)  and  ;=J(sr). 

In  order  to  compute  A(s),  we  need  to  determine 
the  distance  by  which  l(s)  has  moved  due  Lo  the 
change  of  contact  at  u  The  cases  when  u>i(s), 
l(s)>u.  or  u  =i(s)  are  updated  separately  as  follows: 

A(s)  =  /„.,(.)(<5(u))  where  u>l(s). 

A(s)  =  /,f,'),u('5("))  where  L(s)>m. 

A(s)  =  6{u)  where  u=i(s). 

These  functions  can  be  computed  by  appropri- 
ate composition  of  the  functions  stored  at  the 
nodes  on  the  path  from  u  to  r.  Functions  at  the 
right  sibling  nodes  along  this  path  are  composed  Lo 
compute  /,j  if  t<u,  and  those  at  the  left  sibling 
nodes  if  i>u.  For  instance,  let  u=7  in  Fig  5.  The 
function  Z,.,  that  is  needed  at  the  node  0  can  be 
computed  by  composing  /to  and  /,  ,  This  estab- 
lishes that  the  function  tree  can  be  appropriately 


updated  in  O(togk)  Lime  afLer  each  conlacl  change. 
Since  Lhere  can  be  0(ti)  conLacl  changes  during  the 
rotalion  of  the  supporting  polygon,  the  total  time 
complexity  of  this  algorithm  is  0(n log*). 

5.    Conclusions 

We  have  considered  the  problem  of  finding  a 
minimal  vertex  polygon  nested  between  two  convex 
polygons,  and  have  oflered  three  algorithms,  which 
find: 

(1)  a  polygon  with  at  most  one  vertex  more 
than  the  mirumum  in  0(71)  Lime; 

(2)  a  minimal  polygon  In  0(71*)  time; 

(3)  a  minimal  polygon  in  0(7ilogfc)  Lime. 

We  suspect  that  the  last  algorithm  is  optimal, 
but  have  not  proven  this.  It  Is  worth  noting  Lhat  the 
technique  of  detecting  contact  changes  and  com- 
puting the  minimum  between  changes  has  been 
successfully  applied  to  several  other  geometric 
extremal  problems  |DA],  [OR],  |To]. 

A  natural  extension  of  this  problem  is  to 
remove  the  convexity  assumption  find  a  minimal 
vertex  polygon  nested  between  Lwo  simple  polygons 
An  0(n^)  algoriLhm  has  reccnlly  been  prescnLed 
[SOJ.  Finally,  Lhe  original  problem  as  posed  by 
Klee.  tor  polyLopes  of  arbiLrary  dimension,  remains 
open  even  in  the  case  of  Ihree  dimensions. 
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Fig.  5      Function  tree. 

(path  of  updates  shown  by  arrowheads  when  u=7) 
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